Combinatorial formula for Macdonald 



o 
o 
o 

aoj polynomials, Bet he Ansatz, and generic 

< ; Macdonald polynomials 

m 

Andrei Okounkov* 

< 

J2 ' Abstract 

We give a direct proof of the combinatorial formula for interpo- 
lation Macdonald polynomials by introducing certain polynomials, 
which we call generic Macdonald polynomials, which depend on d 
additional parameters and specialize to all Macdonald polynomials of 
degree d. The form of these generic polynomials is that of a Bethe 
\ eigenfunction and they imitate, on a more elementary level, the R- 

' matrix construction of quantum immanants 

00 ' 
O 

1 Introduction 

1.1 Interpolation Macdonald polynomials and combi- 
natorial formula 



x 



1.1.1 The polynomials J M 

Consider the symmetric polynomials 

...,x n ; q,t, s) G Q(q,t, s)[xi, . . . ,x n } s{n) , 

where \i is a partition with at most n parts, satisfying the following interpo- 
lation condition: for any partition A such that fi <£_ A 

J M (A) =0, \=([l,\ 1 ],...,[n,\ n ]), (1.1) 
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where 

sqi +1 t i+1 

lhj\-—^r + — ri - 

Together with the condition that the degree of J M is at most this de- 
termines 1^ uniquely up to a constant factor. These polynomials / M were 
introduced and studied in [jID| ; they specialize to interpolation polynomials 



studied by Knop, Olshanski, Sahi, and the author in a series of papers, see 
the References. We will call the polynomials 1^ the interpolation Macdonald 
polynomials. 

1.1.2 Combinatorial formula 

The polynomials 1^ are very distinguished multivariate special functions with 
a large number of deep properties and important applications (in particular, 
to the better known orthogonal Macdonald polynomials). Most of their prop- 
erties can be quite easily deduced from the following explicit formula for these 
polynomials which is known as the combinatorial formula. 
Consider the diagram of the partition fi and a function 

/:D^{l,...,n} 

from the squares of this diagram to integers. The function / is called a reverse 
tableau if its values weakly decrease along the rows and strictly decrease along 
the columns of \i. We have the following formula (1^ 



/ sg col(D) j/(D)+row(D)\ 

J M= */ 11 ( ~ t /(D)+row(D) - s(f o\{U) ) > ( L2 ) 



reverse tableux / Odfi 



where col(D) and row(D) denote the the column and row numbers a square 
□ G fi and is a certain weight of the tableau / which will be specified 
below. It is the same weight as appears in the combinatorial formula for 
the ordinary Macdonald polynomials, which makes it clear that the ordinary 
Macdonald polynomial is the top degree term of 1^. 

1.1.3 Direct proof of the combinatorial formula 

One of the main results of this note is a proof of the combinatorial formula 



( |1.2j ) which does not use any properties of 1^ except for their definition. We 



recall that a great deal of nontrivial properties of 7^ were used in JED] in the 
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proof of AO]). Now many of these properties can now be without difficulty 
deduced from ( [L.2|) without falling into a vicious circle. 

It is obvious that the right-hand side of ( |1.2| ) is a polynomial of degree 
|ju| and it is also easy to check that it satisfies the vanishing condition ( |1 . 1| ) , 



see for example Lemma 4.1 in [LTj|. What remains is the nontrivial task of 



checking that this polynomial is symmetric in xi, . . . , x r , 
1.1.4 How we prove the symmetry 



Our way of verifying the symmetry in (|1 . 2|) will be the following. We intro- 
duce certain polynomials 

Iu{x) e Q(ui,...,u d , q,t, s)[xi,...,x n ] 

of degree d in x. This polynomials have the property that for any /i such 
that = d one can tune the parameters {ui} so that to make Iu go into the 



right-hand side of (|1.2|) . Because of this property, we call the polynomials 
In the generic Macdonald polynomials. 

The additional freedom which comes with the parameters makes it 
easy to check that the polynomials Iu(%) are symmetric. This is done in Sec- 
tion ||. After that, in Section |3] we show how to specialize the polynomials l\j 
to polynomials in the right-hand side of ( |1.2| ) which, in particular, completes 
our new proof of (|1.2|). 



1.2 Generic Macdonald polynomials 

1.2.1 Shifted Schur functions and Quantum Immanants 

The polynomials l\j are not only a convenient technical tool but also an 
interesting object of study in their own right. 

In order to better motivate their definition, we begin with the special case 
when q — t — > 1 and s — > oo. In this limit, the polynomials J M become (after 
a change of variables) the shifted Schur functions s*, see |T2" |, for which the 
combinatorial formula is 



reverse tableux / Dg/x 



where 

c(D) = col(D) - row(D) 
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denotes the content of a square □ 6 fi . 

The functions s* have a very nice and important representation-theoretic 
interpretations as the image under the Harish- Chandra homomorphism of 
a distinguished linear basis of the center of the universal enveloping alge- 
bra lX(0t n ). The corresponding central elements are known as the quantum 
immanants and denoted by S^. 

Explicit formulas for quantum immanants were obtained in ||, see also 
|5|, [7] . These formulas can be viewed as a noncommutative analog, or quan- 



tization, of the combinatorial formula ( |1.3|) , see for example Section 3.7 of 
||. We recall from H the following construction of S M which was motivated 
by representation theory of the Yangian Y(gl n ). 

Consider the algebra of matrices with entries in U(g[ n ). Let E = (Eij) 
be the matrix formed by the standard basis elements of g[ n . Denote by d 
and consider the following expression 

tr (E - m) ® (E - u 2 ) ® ■ ■ • ® (E - u d ) ■ Yl Rij{ui - %) , (1.4) 

where 

R 1J (n) = l-M G Q( u )^) 



is the rational i?-matrix, the ordered product in ( |1.4| ) ranges over all pairs 
1 < i < j < d ordered lexicographically, the symmetric group S(d) acts on 
the o?-fold tensor product by the usual permutation of indices, and the trace 
means, as usual, the sum of the diagonal matrix elements. 
The expression (|1.4| ) is an element of 

U(Qi n )®Q( Ul ,...,u d ). 



It is explained in how ( |1.4| ) turns into E>^ when the parameters {u^} ap- 
proach, in a certain special way, the contents of the diagram /i. In this limit, 
the product of i?-matrices in ( |1.4j ) turns into a projection onto a certain 
specific vector in the irreducible representation of S(d) corresponding to fi. 



1.2.2 Formula for generic Macdonald polynomials 

The definition of the polynomials l\j is our attempt to keep as much as 
possible from ( |1.4j) in the case of general q and t. This general case has two 
main features. On the one hand, in the absence of a simple structure to 
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replace U(gl n ) for the general q and t, we do not go over the Harish-Chandra 
bridge and stay entirely within the commutative world of polynomials. On 
the other hand, the additional flexibility provided by parameters q and t 
allows us to work with polynomial in a way very similar to what was done 
in the noncommutative world of U(f([ n ). 

As a commutative replacement for the i?-matrix R(u), we introduce the 
following function 

= (1 r qu) [) ~ tu { q) 

HK ' {l-u)(l-tu) 
and define the polynomial Ijj to be the following analog of ( |1.4|) 

n d , ^ t ik \ 
I u (x 1 ,...,x n ) = II ( Xi * + TiT + — ) II P( u k/ui) , (1.5) 

ll,...,ld=l fe=l x 7 lfc<i; 

where the summation replaces the trace, first product replaces the tensor 
product, and the second product replaces the product of R- matrices. 

Of course, the general structure of the formula ( |1.5| ) is a typical form of 
a Bethe eigenfunction. 

1.2.3 Properties of generic Macdonald polynomials 

As already mentioned, in Theorem [TJ below we prove that polynomials l\j are 
symmetric in x±, . . . , x n . The proof of the symmetry is elementary and based 
on the study of the singularities of Ijj in parameters {ui}. 

In Section |3| we show that as the parameters {u^ approach the (q,i)- 
analogs of the contents of a diagram /i, that is, as 

{u t } - {const g coi(n) r row(n)| n ^ (16) 

the generic Macdonald polynomial becomes the interpolation Macdonald 
polynomial labeled by /i, see Theorem [2| The limit ( |1.6| ) has to be taken 
in a certain special way which is parallel to what one does to obtain as a 
limit of flOp . 



Finally, we mention an important but potentially difficult open problem 
which is to expand the generic polynomial ly in terms of the polynomials 1^. 
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2 The polynomials Ijj(X) 



2.1 Definition 

Let X = {xi, . . . , x n } and U = {ui, . . . , Ud} be two finite sets of variables. 
Define a polynomial Iu{x) of degree d = \U\ in the variables X with coeffi- 
cients in rational functions in U as follows. 
Introduce the following function 

(1 -qv)(l -tv/q) 
P{V) = (l-v)(l-tv) ' 

where q and t are parameters. Given a function / 

/:{l,...,d}->{l,...,n} 

set, by definition, 



> Uj ) , 



n p ( ui i i 

where the product is over all pairs i,j in {1, . . .,d} such that f(i) < f(j). 
The polynomial Iu{x) is defined by the following formula 



/ i=l \ 4 / 



(2.1) 



where the summation is over all n d possible functions /. It is clear that 
the definition ( |2.1| ) is just another way of writing the formula (|1.5|) from the 
Introduction. We call these polynomials Ijj the generic Macdonald polyno- 
mials. 

From definition of Iu{X), the following decomposition property of this 
polynomial is obvious. Fix any k — 1, . . . , n and set X' = {x\, . . . ,Xk}- Then 

Iu(X) = J u'(X') I UW (X \ X') J] p( Ui / Uj ) , (2.2) 
u'cu UieU',Ujeu\U' 

where the summation is over all subsets U' C U. 
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2.2 Symmetry 



Since no ordering on U was used, the polynomial Iu{X) is obviously invariant 
under permutations of the Wj's. Our main result about the polynomials Iu{X) 
is that it is also invariant under permutations of the Xj's. 

Theorem 1. The polynomial Iu{X) is symmetric in xi, . . . , x n . 

In the proof, we shall need some elementary properties of the function 
p(v) which we collect in the following 

Lemma 1. We have 

P(v) =p(^)> ( 2 - 3 ) 

p(v) = l+ ( q ~ 1) ( q ~ t) v + 0(v 2 ) , v^O, (2.4) 
q 

^) = l+ te ~ 1) ''~" - + o(-L), „^oo. (2.5) 
qt v \ v J 

Proof of Theorem^. It suffices to check that Iu{X) is invariant under a 
transposition of two adjacent a^'s. By the decomposition property (|2.2|), 
this reduces to showing that Iu(xi, ^2) is symmetric in x\ and x 2 , identically 
in , Ud for any d. 

Let us single out one of the u^s and write 

g(ui) = Iu(xi, x 2 ) - Iu(x 2 , x x ) , 

where x±, x 2 , and the rest of the u^s are considered as parameters. Our goal 
is to show that g = identically. 

The function g is a rational function of u\ with at most first order poles 
at the following points 

U\ E {0, 00, Ui, t ±l Ui] , i = 2, . . . , d . 

First we show that, in fact, this function is regular everywhere. By symmetry, 
it suffices to consider the points {0, 00, u 2 , t ±1 M 2 }. 

Observe that on the divisor u\ = u 2 the function Iu{xx, x 2 ) itself is regu- 
lar. Indeed, a function which is symmetric in u\ and u 2 cannot have a pole 
of exact order one on the diagonal U\ = u 2 . 
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Next, consider the divisor u\ = t 1 u 2 . Only those functions / which 
satisfy 

/(1) = 1, /(2) = 2 

contribute to the pole of Ijj{xi,x 2 ) at u\ = t~ x u 2 . For such function / and 
any i = 3, . . . , d we have either f(i) = /(l) < /(2) or /(l) < /(2) = /(i). 
The condition u\ = t~ l u 2 and fl2.3Q yield that 



«i \ /Ml 
P\ — )=P[ — 



Therefore 



R f ~ R f\{Z,...,d] 



Up 

i=2 



where /|{3,...,d} denotes the restriction of / to {3, ... , d}. It follows that the 
residue of Iu(xi, ^2) at Wi = t~ l u 2 is proportional to 



Ui t \ ( u 2 t 2 \ 

+ T + -) {x* + ¥ + -)l{u 3 ,.,u d} (^) 



which by the condition u\ = t~ l u 2 and induction on the number of UiS is 
symmetric in X\ and x 2 . 

In other words, g{ii\) is regular at U\ = t~ 1 u 2 . By the symmetry between 
Ui and u 2 , it is also regular at u\ = tu 2 . 

Now consider the point u\ — 0. Since for any function / we either have 
/(l) = lor/(l)=2, it follows from the definition of Iu{xi, x 2 ) that 



Iu{xi,x 2 ) = ^ RfUf 
f 



7 /'(0=2 



7 /'(i)=l V 



, (2.6) 



where the summation is over all possible function 

/' :{2,...,d}-{l,2}, 
and where we use the following abbreviation 



% = 11 ^ W + ^ + — j . 



From (373) we have 



n » - 

I Mi/ - LJ - \ Ui 

17 /'(i)=2 V 



t t(q-l)(q-t) 1 ^. . 

— +xi + ^ ^ ^ > - + OM, ui 

Mi O Mj 

H /'(i)=2 

Similarly, from (|2.5|) we obtain 



/'«=1 

^2 



t 2 t(q-l)(q-t) 1 _ . 



til O z — ' Ui 

H /'(*)=! 



Therefore, 

'^ + ^ 2 *(?- NT* 1 \ r / \ /~i / \ 
h Xi + X 2 H 2^ _ J {«2,..,« d } 0^1, «2j + 0{Ui) 



mi g ~ Ui 



as Mi — > 0, which means that both the Mj -1 and u\ terms are symmetric in 
X\ and %2- In other words, not only is g{u\) regular at Mi = but it also 
vanishes at Ui = 0. Similarly, g{u\) is regular at oo and thus, as a function 
regular everywhere and vanishing at one point, g is zero identically. □ 



3 Specialization to interpolation Macdonald 
polynomials 

3.1 Limit transition 

Let A be the diagram of a partition of d. The diagram A has d squares and 
it will be convenient to consider the set U to be indexed by the squares of A 
rather than integers {1, . . . , d}. Denote by I\(X; q, t, s, r) the image of Iu{X) 
under the following specialization of the variables U 

UU = -S g col (n) r r°w(D) ; □ £ A . (3.1) 
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Here s and r are parameters, □ G A is a square in the diagram A, and col(D) 
and row(D) are the column and row numbers of □, respectively. Also, the 
functions / will now be functions from the diagram A to {1, . . . ,n}. Such 
functions are known in combinatorics as tableaux on A. 

Let Rf(X;q,t,r) be the image of Rf under the specialization Q3.1]) . Be- 
cause only the ratios of the tt's enter the formula for Rf, this specialization 
does not depend on s. This image is nontrivial only for some special /'s as 
the following proposition shows: 

Proposition 1. We have 

R f (\;q,t,r) = 

unless f is row decreasing, that is, unless the values of f weakly decrease 
inside the rows. 

Proof. Suppose / is not row decreasing. Then there exist two squares □ and 
□ lcft in A such that n ieft is just left of □ and 

/(□ left ) < /(□) . 

In this case, Rf(X; q, t, r) has the factor of 

^□icft \ ( 1 



and hence -R/(A; q, t, r) = 0. □ 

The main result of this Section is the following: 
Theorem 2. We have 

h{X; q, t, s) cx lim (1 - rt)^- 1 / A (X; q, t, s, r) , (3.2) 

where I\(X; q,t, s) is the interpolation Macdonald polynomial corresponding 
to the partition A and the proportionality constant lies in Q(q,t). 



Theorem |2| will be established by comparing the limit in ( |3.2|) with the 
right-hand side of ( |1.2j ), which in view of the symmetry of the Ijj established 



above, simultaneously proves (|3.2j ) and the combinatorial formula. 
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First, we check that only those / which are reverse tableaux contribute 
to the limit (|3.2| ). Then, in the next subsection, we shall check that for such 
functions / the weights 

R f (X; q, t) = lim (1 - rt) e ^~ 1 Rf(X] q, t, r) (3.3) 

are proportional to the weights \frf in the combinatorial formula (|T 
Proposition 2. We have 

R f (\;q,t) = 

unless f is a reverse tableaux on X. If f is a reverse tableaux then Rf(X] q,t) 
is finite and nonzero. 

Proof of Proposition^. Recall that p(v) has two poles v = 1, i -1 and two 
zeros v = q' 1 , qt' 1 . The pole v = 1 is not an issue in limit ( p.3[) and the zero 
v = q" 1 has been already accounted for in Proposition [I]. The pole v = t^ 1 
means that 

( uu \ 1 

P -»• oo, r -> -, 

\uo u J t 

where D up is the square right on top of the square □. This factor is present 
in Rf(X; q, t, r) if 

/(□) < /(D up ) . 
Similarly, the zero v = qt~ l means that 

P -»• 0, r -> -, 

V M n U p&icft/ 1 

where «□ &left is the square to the left of D up . This factor enters Rf(X; q, t, r) 
if 

/(□) < /(EUwO . 

Let ord/ be the order of the pole of Rf(X; q, t, r) at r = t~ x . We will show 
that 

oid f < i(X) - 1 (3.4) 
and that this maximum is reached precisely for reverse tableaux. We have 
ord; = #{□, /(□) < f(D up )} - #{□, /(□) < /(□„„*«)} . (3.5) 
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The contribution of any square □ G A to ( |3.5| ) is or ±1. 

If the square □ is not in the first column of A then its contribution to 



( |3.5| ) is < 0. Indeed, if /(□) < /(D U p) then since / is row decreasing we also 
have 

/(□) < /(D up ) < /(□ up&left ) ■ 

That is, only the squares in the first column can make positive contribution 
to ( |3.5| ) which proves ( |3.4| ). It is also clear that if / is a reverse tableaux then 
ord/ = £(X) - 1. 

Now suppose that / is not a reverse tableaux. In order to prove that 

ord/ < £(X) - 1 

it suffices to find one square whose contribution to ord/ is negative. By as- 
sumption, there exists a square □ such that /(□) > /(□„„). Additionally, we 
can assume that such a □ is minimal in the sense that f(n ie{t ) < f(O upMeft ). 
Then since / is row decreasing we conclude that /(□) < /(□i C f t ) < /(n up &ioft) 
which means that the contribution of this square □ to ( p.5|) equals —1. It 
follows that ord/ < £(X) — 1. This concludes the proof. □ 

3.2 Computation of Rf( A; q, t) 
3.2.1 The map L 

The weights -R/(A; q, t) are products of a large number of factors of the form 
(1 — q k t l ) where k, I G Z. We will treat such factors as single indivisible 
objects, ignoring the fact that (1 — q k t l ) is reducible if gcd(fc, I) > 1. 

Therefore, it will be convenient to switch, following Macdonald, from the 
multiplicative to the additive notation as follows. Introduce a map L from 
the free Abelian group (written multiplicatively) generated by the symbols 
(1 — a) to the free Abelian group (written additively) generated by symbols 
a which is defined on generators by 

L((l-a)) = a. 

The map L is also a very convenient way to deal with zeros and poles. If 
in an expression A(z) depending on a parameter z 



A(z)=H(l-a i (z)) 



±i 
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m + parentheses in the numerator and m_ parenthesis in the denominator 
vanish at a certain point z = z then L(A(z)) will have a summand of m + — m_ 
at z — zq. 

We agree that 

- <f X ) 

i=0 



-1 



l-q 



which converges as a formal series in q or else converges in the usual analytic 
sense if |g| < 1. Normally, we will only encounter such linear combinations 
of the fractions yz^ which are, in fact, polynomials. 

Taking into consideration the relation (1 — a) = —a(l — a -1 ) we see that 
it is permissible to make the following transformation 

^ia^^ia" 1 (3.6) 

provided n( — a «) ±1 = 1> where the choice of signs is the same as in (|3.6| ). 
Finally, observe that 

LWl))=I (£zi&Z«>, (3.7) 

Upts-D-x ''- 1 ^-" . (3.8) 

These two formulas are related by the transformation ( |3.6| ). 

3.2.2 Combinatorial formula for interpolation Macdonald polyno- 
mials 



The combinatorial formula for the interpolation Macdonald polynomials flO 
expresses the polynomial corresponding to the partition A as a sum over all 
tableaux / on A as follows 



11 ( X /( Q ) ~ t /(D)+row(D) Sf7 col(D) J 



(3.9) 



reverse tableux / DsA 



Here \&/ are a certain weight of associated to / which will be defined mo- 
mentarily. Note that ( |3.9| ) differs by a change of variables from the formula 
established in [[[(] . 
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Recall that a reverse tableaux on A with entries in {l,...,n} can be 
encoded as a sequence of diagrams 

A = A D Ax d A 2 D ■ ■• D A n = 0, (3.10) 

such that each skew shape \k/\k+i is a horizontal strip. This encoding is 
simply the following: 

X k = f-\{k + l,...,n}). 

By definition, 

n-l 

*/ = n ^^Mfe+i ' 

fc=0 

where (see Example 2(b) in Section VI. 7 of ||) 

i<i<j<t(n) 

One can check that this is in fact a polynomial in q and t. 
3.3 The weights R f (X;q,t) 

It is clear from their construction that the weights Rf(X; q, t) factor similarly 
to the factorization of 

n-l 



Rf(X;g,t) = Y[p\ k /\ k+1 , 



k=0 

where 



l (Pa/«) = 

(q-l)(q-t) 



? col(H)-col(D) f ow(D)-rowP) + £ ( A ) _ ^ _ (3 n ) 

□eu, HeA/u 



Here the number £(X) — £(yu) e {0, 1} is precisely the compensation for the 
singularity of -R/(A; q, t, r) as r — > t~ x . 

Proposition 3. The ratio -R/(A; q,t)/^ / depends only on the diagram A and 
not on t/ie particular choice of a tableau f on A. 
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It is clear that this proposition is all what remains to prove in order to 
complete proof of Theorem 

Proof. Let us split the sum in ( p.ll|) in two 

L(Pa//J = £+ + £-, 

where S + and £_ contains summands with nonnegative and negative number 
row(D) — row(KI), respectively. 
Using the identity 

E E * k - 1 = -tt^ (<r c + * b - d - « a - d ~ x - ^ c+1 ) 

k=a l=c \ q > 

to perform the summation along the rows, we compute: 

S + = V (q x '-^ + - q x ' - q^~^) t j ~* . 

1 — q v ' 

Similarly we compute £_ using the formula (|3jjp instead of ( [3.71) and we 
find that 

— ^ ~ _ (q^i- x i+i _|_ g-H+i _ q~ x i+i _ gW-w+i) _ 

1 — <7 ' ^ 

*<i<^0) 

For any diagram r\ define £, v by 

L(e,) = f^ E + - g™ 1 ) - ^) + 1 . 

i<j<i(v) 

Here the constant 1 — £(r]) compensates for the appearances of q°t° in the 
expansion of j^q~ 71:i+1 P~' 1 which occur precisely when i = j < £(r]). Also set 
by definition 

_ Pa/ m £a 
TX/fl ^A/ M ^ ' 
It follows trivially from this definition that 

^ ~ A 1 1 (3 ' 12) 

•f SA fc=0 
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It is clear that 



l 0a/ m ) 



q-t 

l-q 



Note that by our construction we have £(X) < £(fi) + 1. Therefore the above 
sum is either zero or else the only possible value of j in it is j = £{X). In the 
latter case we have A J+ i = which means that 



q 



- 1 - 1 



i<£(\) 



In other words, 



L(7ty„) = <J t - q - 1 
l-q t-1 



l-q t-1 

£(\) = e(n) + i . 



It now clear that 



n-l 



fc=0 



depends only on the first and last element in the chain ( |3.10|) which precisely 
means that ( |3.12| ) does not depend on the particular choice of the tableau / 
on A. □ 
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